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Abstract 

Minimum energy configurations in celestial mechanics are investi- 
gated. It is shown that this is not a well defined problem for point- 
mass celestial mechanics but well-posed for finite density distributions. 
This naturally leads to a granular mechanics extension of usual celes- 
tial mechanics questions such as relative equilibria and stability. This 
paper specifically studies and finds all relative equilibria and minimum 
energy configurations for N = 1,2,3 and develops hypotheses on the 
relative equilibria and minimum energy configurations for N ^ 1 bod- 
ies. 

1 Introduction 

Celestial Mechanics systems have two fundamental conservation principles 
that enable their deeper analysis: conservation of momentum and conserva- 
tion of (mechanical) energy. Of the two, conservation of momentum provides 
the most constraints on a general system, with three translational symme- 
tries (which can be trivially removed) and three rotational symmetries. If no 
external force acts on the system, these quantities are always conserved in- 
dependent of the internal interactions of the system. Conservation of energy 



instead involves assumptions on both the lack of exogenous forces and on the 
nature of internal interactions within the system. For this reason energy is 
often not conserved for "real" systems that involve internal interactions, such 
as tidal deformations or impacts, even though they may conserve their total 
momentum. Thus mechanical energy generally decays through dissipation 
until the system has found a local or global minimum energy configuration 
that corresponds to its constant level of angular momentum. This observa- 
tion motivates a fundamental question for celestial mechanics: 

What is the minimum energy configuration of a N-body system with a fixed 
level of angular momentum? 

This paper shows that this is an ill-defined question for traditional point- 
mass celestial mechanics systems. If instead the system and problem are 
formulated accounting for finite density distributions this question becomes 
well posed and provides new light on celestial mechanics systems. 

This is a well motivated adjustment as real systems always have a finite 
density and, hence, any particle in a celestial mechanics system has a finite 
size. Such a physically corrected system has been called the "Full iV-Body 
Problem," as inclusion of finite density also necessitates the modeling of the 
rotational motion of the components, which is not needed for consideration 
of point masses. It also necessitates consideration of contact forces as their 
mass centers cannot come arbitrarily close to each other, as at some distance 
they will rest on each other. Thus, introduction of this finite density correc- 
tion allows the minimum energy configurations for an A'"-body system to be 
explicitly defined and computed for a given level of angular momentum. 

Given this perspective, an interesting problem is to track the absolute 
minimum energy configuration of a collection of N particles as the system 
angular momentum increases from or decreases to zero. This is, essentially, 
an investigation of the celestial mechanics of granular systems as a function of 
total angular momentum. This problem has been shown to be relevant to the 
understanding of solar system bodies, especially among asteroids whose size 
is small enough so that when the components rest on each other they have 
insufficient gravitational attraction to overcome material strength, and thus 
retain the physical characteristics of rigid bodies resting on each other. Due 
to the celebrated YORP effect, which has established that non-symmetrically 
shaped bodies subject to solar radiation will change their spin rates over 
time, this question has several practical applications and has been implicated 
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in how small asteroids form binary systems [T^]. The question of stable 
minimum energy configurations for the N = 2 particle Full Body problem 
has been worked out in detail flEl [20] . and has been verified as a viable 
physical model via astronomical observations of asteroids |l2j . 

This paper presents a number of Theorems that motivate and enable our 
study of minimum energy configurations, and then explore some exact results 
for the cases of = 1,2,3. Some hypotheses for the case ^ 1 are also 
given. 

Theorem 1: A sharper version of the Sundman Inequality is derived: < 
2IhT < 2/pT, where H is the system angular momentum, T is the 
system kinetic energy, Ip is the system polar moment of inertia and Ih 
is the system moment of inertia about a fixed direction in space. This 
form of the Sundman Inequality is the appropriate generalization for 
bodies with finite density distributions. 

Theorem 2: The minimum energy function £ is defined: £ = + U < 
E, where U is the total potential energy of the system and E is the 
total system energy. It is shown that £ = E occurs if and only if the 
minimum energy function is stationary with respect to all variations 
for a given level of angular momentum H and represents a relative 
equilibrium for the system (borrowing from a theorem originally proven 
by Smale). 

Theorem 3: It is shown that the point-mass A^-body problem with a spec- 
ified level of angular momentum does not have either global or local 
minimum energy configurations for all A^ > 3 (borrowing from a theo- 
rem by Moeckel). 

Theorem 4: It is shown that the finite-density A^-body problem has a global 
minimum energy configuration for all A^ > 1 and for all levels of angular 
momentum. 

Theorem 5-8: All minimum energy configurations, both local and global, 
are found for all levels of angular momentum for the finite density cases 
N = 1,2,3. 

Hypothesis 1: It is hypothesized that for all A^ the local and global mini- 
mum energy configurations for a finite density system will either con- 
dense into a single aggregate and spin at a uniform rate, or will lie in a 
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two component system with the orbit and all rotations spinning at the 
same uniform rate. 

The current paper proves a number of fundamental results that motivate 
and enable the further exploration of these questions. These results show a 
surprising complexity in the evolution of minimum energy states as a func- 
tion of angular momentum, with distinctly different pathways arising as the 
number of particles in the system increases. 

2 Background 

2.1 Conservation Principles 

A closed mechanical system has two fundamental conservation principles: 
conservation of momentum, both translational and rotational, and conserva- 
tion of energy. The statement of these is simple and can be derived starting 
from the fundamental equations of motion for this system. This paper does 
not detail these results, as they are well known (see the explicit derivation 
of these conservation principles for the full 2-body problem in [17J). 

The result is that the total translational momentum, angular momentum 
H and total energy E = T + U (where T is kinetic energy and U is gravi- 
tational potential energy) is conserved under motion within the system. By 
specifying the barycenter of the system to be at rest and at the origin of the 
system the three components of the translational momentum are removed, 
which allows the fundamental quantities of the system to be stated in a purely 
relative form. Angular momentum of the system is also conserved so long 
as no exogenous forces act - regardless of the internal actions between the 
components. The same is not true of the energy, as its conservation depends 
on both the lack of exogenous forces and on strict specification of internal in- 
teractions. Specifically, all internal interactions must be reversible, meaning 
that no dissipation can occur. For ideal celestial mechanics involving point 
mass interactions this is a reasonable model - however the point-mass celes- 
tial mechanics model is always an approximation due to its infinite density 
distribution at the particle location. The truth is that all physical celestial 
mechanics systems of interest involve dissipative interactions, and hence all 
natural systems will tend to dissipate energy, meaning that total energy is 
no longer conserved. In many cases the timescale for energy dissipation is 
quite slow, allowing for the system to be adequately modeled over extended 
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periods of time as conserving energy. However, as has been well established 
in the study of the solar system, energy dissipation due to tidal distortion 
of bodies and internal stress wave propagation eventually dominates the dy- 
namics and interactions of almost all celestial mechanics systems and forces 
them to evolve to lower energy states, all the while conserving total angular 
momentum |B1 Hj . 

Two simple examples of this can be mentioned. First is that almost all 
minor planets in the solar system are in or close to uniform rotation about 
their maximum moments of inertia - this despite a past filled with mutual 
impacts and gravitational torques from their occasional interactions with 
planets. This is simply due to the fact that uniform rotation about the 
maximum moment of inertia is the minimum energy state for a solo rotating 
body. The actual dissipation occurs due to the interaction of time varying 
accelerations that the body experiences when not in uniform rotation with 
the material properties of the body, causing minute elastic deformations that 
turn mechanical energy into heat, which is then radiated away [2]. The other 
example is the mutual orbit of the Earth and moon. The moon has already 
shed its excess rotational energy and has settled into a local minimum energy 
state, with the same face always pointing at the Earth. The Earth still spins 
more rapidly than the mutual orbit, however, raising tides on the Earth 
which dissipate energy and slow the Earth's rotation rate. Conservation of 
angular momentum dictates that the mutual orbit expands - this occurs 
physically by the torque that the Earth's tidal bulge places on the moon. 
The eventual state of the Earth- Moon system (ignoring the sun) is for both 
systems to devolve into mutually synchronous orbits, as is the case for the 
Pluto-Charon system. Once in this final state no excess energy dissipation 
can occur as the system is in its minimum energy state for its given total 
angular momentum. 

2.2 Density Distributions and Fundamental Quanti- 
tites 

To discuss minimum energy configurations for the point mass or finite den- 
sity N body problem requires the definition of a few basic scalar and vector 
functions and quantities. Consider the total kinetic energy, gravitational 
potential energy, polar moment, and angular momentum of an arbitrary col- 
lection of mass distributions, denoted as Bi,i = 1,2, . . . , N . Each body 
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Bi is defined by a differential mass distribution drrii that can fall into one of 
two forms, a point mass distribution or a finite density distribution, denoted 

as 



drrii — 



rriiSdr — Tj|) Point Mass Density Distribution 
Pi{r) dV Finite Density Distribution 



(1) 



where rrii is the total mass of body Bi, S{—) is the Dirac delta function, | — | 
denotes the Euclidian norm, pi is the density of body Bi (possibly constant) 
and dV represents the differential volume of the body. If Bi is described by 
a point mass density distribution, the body itself is just defined as a single 
point Ti- Instead, if the body is defined as a finite density distribution, Bi is 
defined as a compact set in over which Pi{r) ^ 0. In either case the Bi 
are defined as compact sets. 

Assume that each differential mass element dmi{r) has a specified posi- 
tion and an associated velocity. For components within a given body Bi a 
rigid body assumption is made so that the entire body can be defined by the 
position and velocity of its center of mass, its attitude, and its angular ve- 
locity. Finally, assume that these positions and velocities are defined relative 
to the system bary center, which is chosen as the origin, or 



Given these definitions an integral form of the kinetic energy, polar mo- 
ment, gravitational potential energy, and angular momentum vector can be 














(3) 



6 



stated as 



T 



I ir-r)dmi{r) (4) 

N 

h = ir.r)dm.{r) (5) 

u = -gj] V / / ^^^^ (6) 

N 

if = ^ / {r X r) dm,{r) (7) 



N 

where Vij = Vj — rj. Note that the definition of U in Eqn. |6] ehminates 
the self-potentials of these bodies from consideration. As the finite density 
mass distributions are rigid bodies this elimination is reasonable. This nota- 
tion can be further generalized by defining the single and joint general mass 
differentials 

TV 

dra{r) = dnii{r) (8) 

i=l 

Af-l N 

dm{r) dm'{r') = ^ ^ drai{r) dm'j{r') (9) 

i=l j=i+l 

and the total mass distribution B = {Bi, i = 1,2, . . . , N}. Then the above 
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definitions can be reduced to integrals over B: 



T = \ I {r-r) dm{r) (10) 
Ip ^ I {r-r)dm{r) (11) 



H = I (rxr) dm{r) (13) 
Jb 

I rdm{r) = (14) 
Jb 

rdm{r) = (15) 



B 



2.3 Point-Mass Results 

Assuming that all of the differential densities follow the point mass density 
distribution, the classical results for the A^-body problem are recovered: 



T 



\^m{ri-ri) (16) 

N 

h = ^miri-ri) (17) 

i=l 

N-\ N 

u - -^E E ^ (18) 



i=l j=i+l 

N 



I 



H = ^mi(r, xr-i) (19) 



i=l 



It is convenient to apply Lagrange's Identity to the kinetic energy, polar 
moment of inertia and angular momentum, allowing them to be restated 
in a relative form. Note that the following relationships only hold if the 
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barycenter of the system is at the origin. 

N-l N 



T 



EE 



TfLiTTlj \ ' ^ ij ) 



i=l j=i+l 
N-l N 



H 



2M 

i=l j=i+l 
N-l N 



M 



(20) 
(21) 
(22) 



1=1 j=i+l 



2.4 Full-Body Results 



If finite density distributions are assumed for each body the point-mass re- 
sults must be generahzed to incorporate rotational kinetic energy, rigid body 
moments of inertia, angular velocities and explicit mutual potentials that are 
a function of body attitude [T7]. In the following the zth rigid body's center 
of mass is located by the position and has a velocity f,. In addition to its 
mass rrii, the zth body has an inertia dyadic /«, an angular velocity vector Qi 
and an attitude dyadic that maps its body-fixed vectors into inertial space, 
Ai. The basic quantities are then defined as (generalizing results from [T7]): 



N 



T 



^ ^ [rrii (fi • fi) + rij ■ /j • O 

1 

E 



i=l 

N 

m 

i=l 

N-l N 



{Ti ■ Ti) + -Trace(/i 



U 



H 



^ij (''^ijy -^ij ) 

i=l j=i+l 
N 

^ [nii {vi X Vi) + Ai li rii 



(23) 
(24) 
(25) 
(26) 



1=1 



In the above the inertia dyadics are all specified in a body-fixed frame and 
thus are constant, the Uij are mutual potentials between two different rigid 
bodies i and j and are only a function of their relative position and relative 
attitude, Aj^, equal to AJ ■ Aj, and which transfers a vector from the body 
i frame into the body j frame. 
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The kinetic energy, polar moment of inertia and angular momentum can 
again be stated in relative form between the center of masses, leaving the 
rotational components in their current form. 



N-l N N 



2M ^ '-" 2 

i=l j=i+l i=l 



^ N-l N ^ N 

^ mJ2Y1 ^^^J ^^^^ ■ ^^j) + 9 IVace( Ji) (28) 



i=\ j=i+l i=l 



H = 



^ N-l N N 

U^Y^ ^i^i (^^J- ^ + YA-Ii-^i (29) 



M 

1=1 j=i+l 1=1 



2.4.1 Alternate System Moment of Inertia 

An alternate version of the system moment of inertia can be defined once the 
direction of the total angular momentum is specified. First define the total 
inertia dyadic of the N body system as 

I — — ■ rdm (30) 
Jb 

where r is the cross product dyadic defined such that axb = db = ab. 
When evaluated in detail, and in a common inertial frame, the system inertia 
dyadic is equal to 

N 

I = J2['^d^iU-riri) + A-Ii-Aj] (31) 

i=l 

where the identity —a b — {a ■ b)U — ba is applied and where U is the unity 
dyadic and the multiplication of two vectors is defined as a dyad. Given a 
defined direction of the total angular momentum vector in inertial space, H, 
the moment of inertia relative to this direction can be defined as 

Ih = H I H (32) 

This alternate version of the polar moment of inertia is of use in developing 
sharper limits for the Sundman Inequality later. For the point mass density 
distributions all of the moments of inertia disappear and 1^ — Ih if the bodies 
and their velocities all lie in a common plane. 
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2.4.2 Minimum Distances between Mass Collections 

A crucial aspect of the finite density distributions are that these bodies can 
come into contact with each other without singularities or deformations. 
Thus, under a rigid body assumption they will have minimum distances 
between their centers of mass as a function of their relative attitude. If the 
bodies are convex, strong results on the minimum distances between the bod- 
ies can be found. However, even complex shaped bodies will still have well 
defined minimum distances between the mass collections. 

Each body has a minimum diameter (equal to the diameter of the in- 
scribing sphere and denoted by d) and a maximum diameter (equal to the 
diameter of the circumscribing sphere and denoted by D), these are equal 
only if the body is a sphere. The relative distance between any two mass 
distributions Bi and Bj with minimum diameters di and dj is bounded from 
below by the average of these, or > ^{di + dj) = dij. For distances 
Tjj > ^{Di + Dj) = Dij the relative position and attitude between the two 
collections is unconstrained. For any combination of relative position DijVij 
and relative attitude Aij the distance between the two mass distributions 
can be decreased to their relative distance at which they will touch, defined 
as dij{rij, Aij). The contact points of bodies i and j that define this dis- 
tance lie in their respective boundary sets of the mass distributions. For the 
given definition, the set of these minimum distances is a closed set and is 
bounded from below by dij, the minimum distance. For non-convex shapes 
these minimum distance sets can be complicated and may not even be able 
to attain the theoretical minimum. Also, in these cases different definitions 
of the minimum distance can be defined, depending on whether the bodies 
are star-convex or have a more convoluted shape. Note that across any of 
these shapes it is always possible to define a minimum bound on the dis- 
tance between these shapes that can be realized, meaning that the distance 
between two mass distributions is a closed set along its lower bound. If the 
bodies are convex this follows trivially. 

2.4.3 Finite Density Sphere Restriction 

This paper mainly focuses on the sphere-restriction of the Full-Body problem, 
where all of the bodies have finite, constant densities and spherical shapes 
defined by a diameter di. This allows for considerable simphfication of the 
mutual potentials, although the rotational kinetic energy, moments of inertia 
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and angular momentum of the systems are still tracked. In this case the 

moment of inertia of a constant density sphere is rriidf/ 10 about any axis 
and the minimum distance between two bodies will be dij = {di + dj)/2. The 
resultant quantities for these systems are 

^ = ^ E E (^^.- • + 1 E (33) 

i=l j=i+l 1=1 

N-1 N 

^ = -5E E ^ (34) 

i=i j=i+i I ! 

1 ^ m d^ 

^ = M E E '^i'^i ^ ^'i) + E if"' 

1=1 j=i+l i=l 

with the two different versions of the moment of inertia 

N-l N „ AT 



^ E E (^^j- • ^ij) + ^ E IF ^^^^ 

i=l j=i+l 1=1 

'« = ]g E E '"""i (4 - (» • + E ^ (37) 

i=l j=i+l i=l 



2.5 Energy Dissipation Interaction Models 

Implicit in these discussions, although not explicitly incorporated into the 
interaction models, are the dissipative effects of surface Coulomb friction 
and the tidal distortion of gravitationally attracting bodies. These physical 
effects serve dual purpose, in that they will tend to synchronize collections 
of bodies, either resting on each other or orbiting each other, and will also 
dissipate excess energy in the system. 

No surface interaction models between the finite density bodies is con- 
sidered other than Coulomb friction between surfaces. This is needed in 
order for a resting collection of particles to dissipate relative motion between 
each other when in contact and thus also represents one possible mode of 
energy dissipation. Inclusion of this notional model ensures that contact 
configurations will, when reduced to their minimum energy state, all rotate 
at a common rate. It is possible to explicitly incorporate additional surface 
potentials between particles in close proximity, such as the Lennard- Jones 
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potential that models van der Waals cohesive forces [3]. These extensions 
are not considered here, however. 

It is also possible to dissipate energy and synchronize spin rates even if 
bodies are not in contact. Tidal distortions arising from relative motion be- 
tween gravitationally attracting bodies will also transfer angular momentum 
across the system and cause the dissipation of energy. Even if these effects 
are small, as they can be for asteroidal bodiesfT], they are pervasive and will 
cause continual dissipation of energy for systems that are not in a relative 
equilibrium state. 

The ubiquity and pervasiveness of energy dissipation in the solar system 
and its role in the long-term evolution of bodies of all sizes motivates the 
main question concerning the minimum energy states for celestial mechanics 
systems at a given value of angular momentum. 

3 The Sundman Inequality, Amended Poten- 
tial, and Minimum Energy Function 

3.1 The Sundman Inequality 

The fundamental result for this study is found by application of the Sundman 
Inequality. For generality it is applied to the integral form of the angular 



momentum vector, given in Eqn. 13 As this is not the usual form of the 



equations, a descriptive proof is given that relies on the Cauchy-Schwarz 
Inequality. 

Theorem 1. < 21 hT < 2/pT across all of the differential mass formu- 
lations defined in Eqn. [I] and system moments of inertia defined in Eqns. \Tl 
andlM 

The outermost inequality is the usual Sundman Inequality, but the sharper 
limits are new and are distinct for the full body problem. 

Proof. First prove the inequality Ih ^ Ip, which establishes the ordering on 
the right. First ignore the center of mass terms and only consider the mo- 
ments of inertia. Taking each term independently leaves H ■ Ai- I-i- AJ ■ H < 
|Trace(/j). If /i < /2 < I3 are the principal moments of inertia of /j, then 
H ■ Ai - li - AJ ■ H < I3. Thus the inequality reduces to < ^ (Ji + I2 + I3) 
or Is < Ji + /2. However, this inequality and all of its permutations are a 
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fundamental property of mass distributions and moments of inertia. Next, 
consider the center of mass terms, leading to the following inequality for each 
term: — {H ■ < r^, which can be trivially shown to be true. 

To finish, consider the inequality < 2TIh, which is not the usual 
Sundman Inequality, using the general mass integral form of the angular 
momentum. 

First recall that H = HH = J^r x r dm. Dotting both sides by the 
(constant) unit vector aligned with the angular momentum vector yields the 
equality 



H = 



H ■ {r X r) dm 



r ■ (H X r) dm 



Now apply the triangle inequality to the integral to find: 



r ■ {H X r) dm < 



fin 


H X r 


Jb 





dm 



(38) 
(39) 

(40) 



Squaring the original term, equal to H^, and applying the Cauchy-Schwarz 
inequality yields the main result. 

1 2 



B 



H X r 



dm 









2 




H X r 


dm 


< 


Jb 









[H X r) ■ {H X r) dm 



dm 



(41) 



(42) 



But /g dm = 2T and /^(JT x r) ■ {H x r) dm ^ H ■ j^-r ■ r dm ■ H ^ 
H I H^Ih Thus, < 2TIh. □ 

Note that for the point mass density distribution Ip — Ih if all of the 
bodies lie in a single plane perpendicular to H, as the moments of inertia 
are identically zero for a point mass. 



3.2 The Minimum Energy Function and the Amended 
Potential 

The Sundman Inequality provides an important, and sharp, lower bound on 
the system energy for a given angular momentum. The derivation of this is 
simple, but the result has not been extensively used. 



14 



Theorem 2. The total system energy E = T + U is bounded below by the 
minimum energy function, defined as £ = + U , which is only a function 
of the system total angular momentum and the relative configuration of the 
components within in the system. Thus, given a total angular momentum 
for the system, H , the system energy is constrained by £ < E. Equality of 
the system energy and the minimum energy function occurs if and only if the 
minimum energy function is stationary with respect to all possible variations, 
and corresponds with the system being in a relative equilibrium. 

Proof. In the updated Sundman Inequality, < 2TIh, replace the kinetic 
energy with T = E — U , where E is the total system energy and U is the 
gravitational potential of the system. Rearrangement of the terms yields the 
result + U < E, and also defines the minimum energy function. From 



Eqns. 32 and 25 note that the minimum energy function is only a function of 
the system configuration relative to itself, and does not involve any externally 
defined reference points. 

The relative equilibrium result arises from a theorem by Smale |22l |23] , 
proven in a more direct manner by Arnold in [1]. First, note that the min- 
imum energy function is related to the "amended potential" (with possible 
differences depending on the form of the polar moment of inertia used) which 
has found many uses in the derivation and development of constraints on 
motion in the point mass A^-body problem. In [23] it is proven that the 
stationary points of the amended potential correspond to relative equilibria 
of the corresponding point-mass A^-body system at a given level of angular 
momentum. In [21] the proof is extended to more general dynamical systems 
with symmetry, which contains the Full Body system (as established in 

To finish, equality of E and £ is proven at a relative equilibrium. A 
collection of bodies in a relative equilibrium will all rotate at a constant rate 
and maintain a constant polar moment of inertia [H]. This implies that the 
angular momentum vector is an eigenvector of the locked inertia matrix and 
the value In is an eigenvalue of this matrix (see also [Ml [lOl HE])- Thus the 
total angular momentum is H = I ■ fl = Ih^H and the kinetic energy is 
T = ^fl ■ I ■ fl = ^Ih^^- Substituting into the minimum energy function 



"'^Note that the initial, direct appUcation of the Cauchy-Schwarz inequaUty yielded a 
polar moment of inertial Ip. If this moment were used to derive the minimum energy 
function it would no longer equal the amended potential except for the special case of a 
point-mass system. The modified version of the polar moment, Ih is seen to precisely 
yield the amended potential when finite density considerations are taken into account. 
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yields E = PQ.'^/{2Ih) + + U ^T + U ^ E. 

Now consider the case when the system is not in a relative equilibrium yet 

TJ-2 

the minimum energy function equals the system energy, ^j— + U = E. Since 
the system is not in a relative equilibrium, the minimum energy function 
(amended potential) is not at a stationary value. Thus, it is possible to 
take an allowable variation in the configuration to increase its value, or £ + 
S£ > £ = E. However, then the Sundman Inequality is violated, forming a 
contradiction. □ 

CoroUciry 1. The energy is bounded by an additional function as: ^ + U < 
^ + U<E. 

Proof. This follows immediately from the inequality Ih ^ Ip- n 

3.3 Existence of Global Energy Extrema 

Given the defined minimum energy function a rigorous approach can be for- 
mulated to the original question of what the global minimum energy config- 
urations of a celestial mechanics system are. This, it turns out, is strongly 
dependent on the density distribution assumed. The following two theorems 
give explicit results for point mass and finite density distributions. 

Theorem 3. For the Point Mass Density Distribution: 

i) the minimum energy function is undefined for N = 1, 

ii) the minimum energy function has a unique global minimum for N = 2, 
Hi) the minimum energy function does not have a global minimum, or even 
a local minimum, for N >3. 

Proof, i) For A^ = l,Jp = 0, f/ = and the angular momentum is identically 
zero for the assumed barycentric coordinate frame. Hence the minimum 
energy function cannot be constructed. 

ii) For N — 2 a, purely constructive proof is given. At A'^ = 2 the minimum 
energy function terms are 

r - ^^^^-r^ (43) 



mi + 777.2 



U = (44) 
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leading to 



^ _ H^{mi + 7712) _ Qmim2 
2mim2r2 r 

and is defined completely by one degree of freedom, the distance between 
the two mass points. Figure [T] shows a generic graph of this function, and it 
is clear that only a single extremum exists. Taking the variation of £ only 



Generic Energy-Distance Curve 
for a given Angular Momentum 





Single Minimum-Energy 
Configuration at each Angular 
Momentum 



Generic Energy-Angular 
Momentum Curve 



Figure 1: A generic graph of the minimum energy function as a function of 
the distance (left) and of the angular momentum (right). For each angular 
momentum there is a single extrema, a minimum, which corresponds to a 
unique circular orbit. 



involves this term and yields 



5£ 



H'^{mi + 7712) 
7ni7n2r^ 



-1 + 



Q{7ni7n2) r 



6r (46) 



iJ2(mi -|- 7712) 

Setting this to zero for all 6r leads to a unique solution for r* G (0, 00) 

. H'^(7ni + 7712) 

r = 

^(77117712)^ 

The corresponding energy is 



(47) 



^ g^(mim2)^ , 
2H^{7ni + 7n2) ^ ' 

and becomes arbitrarily large as — ?■ and goes to as — )■ 00. Figure 
[1] shows a generic plot of this relation, note that at every angular momentum 



17 



there is only a single energy configuration. To test for a minimum take the 
second variation and substitute the nominal solution to find 



which is strictly positve, and hence the relative equilibrium is at least a local 
minimum of the energy function. To prove that this is a global minimum, 
it can be shown that £ — £* > for all r. To establish this note that this 
inequality defines a quadratic equation in 1/r which can be explicitly factored 
to show that £* is a global minimum at a given value of angular momentum. 
In terms of usual orbital mechanics this global minimum is equal to a circular 
orbit. 

iii) For > 3, note that the relative equilibrium are all stationary values 
of and that £ = E sX the relative equilibrium. By definition these are 
all central configurations with a specific system rotation rate so that they 
remain in a relative equilibrium, and have been analyzed extensively by many 
authors. Moeckel provides a proof that the second variation of the energy at a 
fixed value of angular momentum is indefinite for all central configurations of 
the iV > 3 problem [llj. Thus, there always exist both positive and negative 
eigenvalues of this second variation, meaning that relative equilibrium are 
never even local minima of the minimum energy function, which is equal 
to the energy at a fixed level of angular momentum at a stationary pointj^ 
This result removes any and all central configurations and related relative 
equilibria from consideration as minimum energy configurations. 

To completely establish that there are no global minimum it can be shown 
that the minimum energy function is unbounded from below for any given 

G (0,oo). Again this can be done by construction. For any > 3 
consider a random spatial distribution of the point masses, but allow two of 
them, mi and m2 to be in close proximity, so that their relative distance is 

= e. Then the minimum energy function can be split into two terms, 
£ = £' — Smrm2 ^ ^i^QYB £' IS finite as long as at least one pair of bodies are 
not arbitrarily close to each other. Then, as the bodies are point masses let 
e — to find explicitly that £ ^ —oo. Thus the minimum energy function 
is unbounded from below for an allowable configuration and hence cannot 
have a global minimum across all possible system configurations. □ 

^This is not surprising, as if the second variation were definite it would be trivial to 
constructively prove stability for at least some relative equilibria - implying that the KAM 
theorem would not be needed. 




(49) 



.3 



18 



For a mechanics system, this lack of a global minimum energy configura- 
tion is problematic, as all mechanical systems are expected to have minimum 
energy states. For example, this implies that for the three body problem with 

> there are two final states in the presence of mechanical energy dis- 
sipation. First is that one body is ejected, leaving a two-body system and 
a solo body, with at least one with a well defined minimum energy state. 
Second is that the system remains bound and, under energy dissipation, two 
of the components eventually spiral towards each other until a singularity oc- 
curs. This latter system is unacceptable in a Newtonian system, as the body 
interactions will go beyond the validity of Newtonian gravity and physics at 
some point during their final spiral. 

The problem with this picture lies in the point mass density distribution 
traditionally used for celestial mechanics systems. Now consider minimum 
energy configurations for the finite density distribution. 

Theorem 4. The minimum energy function for an N-body finite density 
distribution has a global minimum energy state for every level of angular 
momentum and for all N >\. 

Proof. This theorem can be proven by showing that all of the mutual po- 
tentials Uij and H'^ / {2Ih) have global extrema across their domains. Then 
the sum of these functions will have global extrema across their combined 
domains. 

For N = 1 mutual gravitational potentials are undefined (as self poten- 
tials are neglected) and C/ = 0. 

For N > 2 mutual potentials between each of the bodies in a collection 
are continuous functions of their domain, which is the relative position vector 
rij and the relative attitude, as can be specified by the transformation dyadic 
Aij. Let us break up the relative position vector into the unit vector and 
its magnitude, Vij — r^jVij. Then the elements fjj and Aij are defined over 
compact spaces, since these live in SO (2) and SO (3). For small values of 
there may be additional constraints on the allowable contact locations and 
orientations between the bodies, but these are all continuously connected 
to the domain of interest (per the earlier definitions) and define closed sets. 
For Ty, the lower limit of separation was defined earlier to be dij. Thus the 
domain r.ij G [dij, oo), and is only closed at its lower end. Also note that all 
mutual potential reach a finite limit of for Vij — )■ oo. Thus, a one-point 



19 



compactification can be made to make the interval over which r^j is defined 
compact. Note that Uij is a well defined and continuous function across this 
domain. Then by the extreme value theorem Uij has a global minimum and 
maximum. 

Now consider the function H^/{2Ih)- The moment of inertia is always 
> for a finite density distribution, as even when = 1 it equals the 
moment of inertia about the direction H, always non-zero by definition. As 
any rij — )■ oo the function 1//// — )■ 0, and thus is well defined. As before, one 
can argue that 1/Ih is continuous over its interval of definition and that its 
domain is compact (or can be made compact). Thus, by the extreme value 
theorem H^/ {2Ih) has a global minimum and maximum. 

To complete the proof, trivially note that the addition of continuous func- 
tions defined over a common compact domain remains continuous. Thus 
the extreme value theorem applies to the minimum energy function S = 
^ + '^f=i^ '^f=i+i ^ij^ it ^ global minimum and maximum. □ 



4 Minimum Energy Configurations in the 
Spherical Full Body Problem 

Having proven existence for the main result, its implications can be explored 
using the simplest possible extension of the point-mass celestial mechanics 
problem to the finite density case, which is the spherical full body problem. 
See previous works by |2H [TOl [T71 UHl El |2Q] for more detailed Celestial Me- 
chanics discussions of the Full 2-Body problem for non-spherical bodies. The 
following analysis strives to be more complete by analyzing all possible op- 
tions, enabled in part by the simpler model of interaction that the spherical 
restriction allows. 



4.1 N = l 

The mutual potential is identically zero and the moment of inertia of the 
sphere is equal through any axis through its center. Thus the minimum 
energy function takes on the simple form 
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and has no degrees of freedom. Thus this system is always in a minimum 
energy state. Had the sphere restriction not been apphed, then the minimum 
energy state would have the body reoriented so that the maximum moment 
of inertia would lie along H |2j . 



4.2 N = 2 

Now revisit the two-body problem for the case of finite density, focusing on 
relative equilibrium, local minimum configurations and globally minimum 
configurations. This problem shows that extensions to finite density mass 
distributions fundamentally changes the energy structure of this problem 
and yields the non-intuitive result that there can be multiple circular orbits 
at a given angular momentum and that there can be unstable circular orbits. 
For generality in the following the two bodies have different masses, mi and 
m2, and diameters, di and 

Theorem 5. There are 0, 1 or 2 orbital relative equilibria for the 2-body 
problem as a function of angular momentum. When there are 2 orbital equi- 
libria, one is energetically stable and the other is energetically unstable. 

Proof. The potential and polar moments of inertia are: 

U = 
Ih = 

Is = 
r > 

Normalize this system by 

r 
H 

Is 

U 



Gmim2 

(51) 

r 

-^r^ + Is 52 

nil + "^2 

^ [rriidl + m2dl) (53) 
= di2 (54) 



r/di2 (55) 
Hy/nii + 1712 



^/Qdi2mim2 
mi + m2 J. 

dl2mim2 ^ 

Udi2 
Qmim2 



(56) 
(57) 
(58) 
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In the following the (— ) notation will be dropped with the assumption that 
all quantities are normalized as given above. With this normalization the 
range of allowable distances between the particles is defined by r > 1. 
The minimum energy function is then 

This is a single degree of freedom function so relative equilibria are found by 
taking a variation with respect to r 



5£ 



1 



r5r (60) 



Setting the variation equal to zero yields an equation for r: 

- + 2Isr'^ + = (61) 

(62) 

Applying the Descartes rule of signs the polynomial either has 2 or no 
solutions, a drastic change from the point-mass problem which always had a 
single solutiorj^ Of particular interest is the point where the transition from 
zero to 2 solutions occurs, which will be a double root of this polynomial. 
Set the derivative to zero to find 

(4r2 - 3H^r + 4/,) r = (63) 

Evaluating both equations simultaneously shows that the bifurcation from 
zero to two solutions occurs at 

r* = VsTs (64) 

Note that the minimum value of Is occurs when mi = m2 which leads to 
di = d2 under an equal density assumption, which leads to Is > 0.4. Thus, 
r* > -\/L2 > 1 for all mass values with constant density and the bifurcation 
structure for = 2 is robust and always occurs when the bodies are separate 
from each other. 



^If — > the single solution from the point-mass case is recovered. 
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An alternative approach is to solve for as a function of r, yielding 



= ' ' (66) 



r 



This has a minimum point, corresponding to the bifurcation of the roots to 



the equation given above in Eqns. 64 and 65 For higher values of H the 



system will then have two roots, one to the right of r* (the outer solution) 
and the other to the left (the inner solution). At the bifurcation point the 
system technically only has one solution. 

The stability of each of these solutions is determined by inspecting the 
second variation of S with respect to r, evaluated at the relative equilibrium. 



Taking the variation and making the substitution for H from Eqn. 66 yields 

{6rf 



5'£ 



(r2 + /,) r2 



(67) 



Checking for when > yields the condition > 3/^. Thus, the outer 
solution will always be stable while the inner solution will always be unstable, 
with the relative equilibria occurring at a minimum and maximum of the 
minimum energy function, respectively. Figure |2] shows characteristic energy 
curves and the locus of equilibria for different levels of angular momentum 
for the spherical full 2-body problem. This should be contrasted with the 
energy function for the point mass 2-body problem, shown in Fig. [T| which 
only has one relative equilibrium. 

□ 

Theorem 6. There are or 1 unique resting relative equilibrium config- 
urations as a function of angular momentum. When it exists, the resting 
equilibrium is energetically stable. 

Proof. There is only a single unique contact configuration for the N = 2 case, 
with the two bodies touching with their centers of mass a distance apart. 
In normalized coordinates the minimum energy function is £^ = n/ r \ — - 
where r = 1 signifies contact between the bodies. For a contact structure 
to be a resting equilibrium, any variations in its degrees of freedom must 
conform to the rigid body constraint (i.e., Sr > in this case) and must only 
cause an increase in the minimum energy function, 6S > 0. If a decrease 
in energy occurs, the system will trend away from the configuration along 



23 




Energy curves at varying 
levels of angular momentum 



1.2 

Separation r 



Figure 2: Locus of orbital equilibria across a range of angular momentum 
values. Plotted is the minimum energy function versus the system configu- 
ration, r, the distance between the two particles. This plot assumes equal 
masses and sizes of the two particles. For clarity, equilibrium solutions are 
shown below the physical limit r > 1. 



an allowable variation in the degree of freedom. Taking the variation in the 
minimum energy function at the contact configuration of r = 1 yields 



5£ 



1 - 



+ 



5r 



(68) 



Since 5r > is the only variation allowed, the relative equilibrium will exist 



when 1 



> 0. Thus if H < 1 + Is the static resting equilibrium will 



exist and ii H > 1 + Is there will be no resting equilibrium. 

When the resting equilibrium exists, any variation in the degree of free- 
dom will only increase the energy, meaning that its variation is positive 
definite and this relative equilibrium is energetically stable. □ 

Theorem 7. The complete bifurcation chart of relative equilibria, minimum 
energy states, and global minimum energy states of the sphere restricted N = 
2 full body problem as a function of angular momentum follows. The sequence 
is graphically represented in Figure\^ 



24 



< < -^^y/h The minimum energy configuration, and the only relative 
equilibrium, is for the two bodies to be resting on each other with r = 
di2- 



v^VTs < H'^ < ^^rr^T"^ Two orbital relative equilibria bifurcate into exis- 
tence at a radius r* = ^/3Is and grow to larger and smaller values of 
r. The outer solution is at a local minimum of the energy, and thus 
is stable. The inner solution is at a local maximum of the energy, and 
thus is unstable. The global minimum energy configuration remains the 
contact solution. 

^^=7^1= < < (1 + The situation is the same as above, except now 
the outer orbital solution becomes the global minimum, and the resting 
solution becomes a local minimum. 

(1 + IsY = The inner orbital solution collides with the resting solution 
and they annihilate each other. This can be called the "fission" condi- 
tion for the resting configuration, as it is the point beyond which the 
resting condition no longer exists and the particles that were resting on 
each other are forced to enter orbit about each other. 

(1 + IsY < Beyond the fission limit there is only one relative equilibrium, 
the outer orbital solution, which is the minimum energy configuration. 










AE = Excess Energy 
• ■ ■ ■ > Increasing AM 
> Decreasing AM 







Figure 3: The energy-angular momentum diagram of relative equilibria (left) 
and the transition diagram for increasing and decreasing angular momentum 
(right) in the Full 2-body problem for equal sized bodies. 
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Proof. Theorems [s] and |6] outline the values of at which the relative equi- 
libria exist and what their stability properties are. The main items to prove 
in the above are i) the transitions between the global minimum configura- 
tion and ii) the collision of the unstable orbital equilibrium with the resting 
equilibrium. 

i) To prove this find the value of that gives equality of the energy of 
the resting and orbital configurations and verify the inequalities. To under- 
stand the global minimum transitions compare the energy of the resting and 
the stable outer orbital configurations with each other. Equating the orbital 
and the resting minimum energy functions at the same level of angular mo- 
mentum, it is possible to solve for the angular momentum at which the two 
functions are equal: 



2(l + /,)(r2 + /, 



H' = ' ^' " (69) 

r(r + l) ^ ' 

To identify the transition point between global minima substitute for the 



angular momentum as a function of orbit radius (Eqn. 66) and simplify the 
resulting equation to find 

(r - 1) P - 2/,r - /,) = (70) 

with an obvious root at r = 1 (discussed in a moment) and a positive real 
root at = Jg + a//s(/s + 1). When r = the global minimum will switch 
from the resting configuration to the outer orbital configuration as angular 
momentum increases. Substituting into Eqn. |66] and simplifying yields 

Hi = i^Oi^L (71) 

Next check the inequalities < ^^^^^^J^ < (1 + hf- First note 

that if Is becomes large the inequahties devolve into 
and strictly hold for Ig > 64/27 for the lower inequality and Is > 2 for the 
upper. Also, taking the minimum value of Is = 0.4 and inserting it into 
the relationships yields 1.9474... < 1.9506... < 1.96. To verify over the 
remaining interval it suffices to plot and compare the functions over this 
compact set. 

ii) Note that the fission of the resting configuration will occur when the 



inner, unstable relative equilibrium touches the value r = 1. From Eqn. 70 
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note that the energy of the resting and orbital relative equilibria are equal at 
this point. Evaluating the angular momentum at this point then yields the 
value H'^ = {l + hf. □ 



Corollary 2. The N = 2 transition diagram as a function of angular mo- 
mentum is discontinuous and exhibits hysteresis. To track the evolution of 
a system under changing angular momentum will require dynamical analysis 
at specific transition points. 

Proof. In Figure |3] the energy of the relative equilibria is graphed when they 
are stable as a function of angular momentum for equal-sized bodies (with 
qualitatively similar results holding across all values). Under increasing or 
decreasing angular momentum the configuration will remain in a locally sta- 
ble relative equilibrium after the global minimum configuration has switched 
to a different configuration. Thus, when the local relative equilibrium be- 
comes unstable or ceases to exist, there is a finite energy difference between 
the current system state and its minimum energy state. If angular momen- 
tum is increased beyond this point this directly implies a period of dynamical 
interaction during which energy dissipation may occur, eventually leading to 
the minimum energy state. 

At fission the system has a corresponding energy of ^(/^ — 1) and thus at 
fission the energy can be positive if > 1. For constant density spheres this 
will occur when the ratio between the masses is < 0.204 ... or > 4.902 . . ., 
assuming constant density [20]. In [17] it is proven that full body systems 
can escape if their free energy (equal to the energy as defined in this paper) 
is positive. Thus, across the range of different mass ratios it is even possible 
for the resulting fission system to eventually escape prior to sufficient energy 
dissipation for being trapped into the outer orbital configuration. In [12] this 
same mass ratio is found in a class of natural bodies that are expected to 
have undergone fission in the past. □ 

For the Spherical Full 2-body problem, these diagrams are notional as the 
actual dynamics of these systems consist of purely elliptic motions. Indeed, 
without explicit energy dissipation the excess energy can never be depleted 
and the system will not evolve. Also, the rotational motion of the spheres 
are completely decoupled from orbital motion in the rigid body limit. How- 
ever, these diagrams present the underlying energy structure for fixed levels 
of angular momentum. Even the slightest departure from rigidity will enable 
energy dissipation to occur if perturbed from either of the relative orbital 
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equilibria. If perturbed from the stable equilibrium the dissipation will allow 

the system to settle into its minimum energy state. If perturbed from the 
unstable equilibrium it will either evolve towards the contact relative equi- 
librium, the stable orbital one or may escape before it reaches this limit if 
the total energy is positive. 



This discussion is restricted to bodies having equal sizes and densities. Thus, 
all particles have a common spherical diameter d and mass m. A restriction 
to the planar problem is also made to simplify the analysis, although non- 
planar configurations can also be analyzed. For the equal mass, point-mass 
celestial mechanics problem, there are only two known relative equilibrium 
configurations, the Euler and Lagrange solutions. When finite densities are 
incorporated there are at least 6 orbital configurations and a total of 9 distinct 
relative equilibria. If non-equal size bodies of similar density are considered, 
there will be additional distinguishable relative equilibria, but this analysis 
is left for future work. 

Given this restriction the polar moment and potential energy take on 
simpler forms. 



where m is the common mass of each body and d the common diameter. 
Then rij > d for all of the relative distances. Now introduce some conve- 
nient normalizations, scaling the moment of inertia by md'^ and scaling the 
potential energy by Qm? /d. Then the minimum energy function is 



4.3 AT = 3 




(72) 



(73) 



ri2 rsi 



£ = 




(74) 



2[(r?2 + i3 + ii)/3 + 0-3] Lfi2 f23 rai 



where 



£ = 



£d 



(75) 




(76) 



Qm?d 
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and the constraint from the finite density assumption becomes fij > 1. 

In the following the — notation is dropped for r^j and H, as it will be 
assumed that all quantities are normalized. 

Theorem 8. There are a total of 9 unique (under symmetry transformations) 
relative equilibria that can exist in the N = 3 sphere restricted full body 
problem. These are shown in Figure^ and described below. 

1. There are two static resting relative equilibrium configurations, the Eu- 
ler Resting and Lagrange Resting configurations. The Lagrange Resting 
configuration is energetically stable whenever it exists while the Euler 
Resting configuration transitions from unstable to stable as a function 
of angular momentum while it exists. 

2. There is one dynamic resting relative equilibrium configurations, the 
"V" configuration, which is always energetically unstable when it exists. 

3. There are 4 mixed orbital and resting relative equilibrium configura- 
tions, the inner and outer aligned and transverse mixed configurations. 
The 2 inner configurations and the outer transverse configurations are 
always energetically unstable. The outer aligned configuration is always 
stable. 

4. There are 2 purely orbital relative equilibrium configurations, the La- 
grange and Euler configurations. These are always energetically unsta- 
ble. 



Proof. 1. Static resting configurations 

First consider the static resting configurations, defined as when all of 
the bodies are in contact and maintain a fixed shape over a range of 
angular momentum values. When in contact there is only one degree of 
relative freedom for the system, defined as the angle between the two 
outer particles as measured relative to the center of the middle particle 
and shown in Fig. |5j As defined the angle must always lie in the limit 
60° < ^ < 300°. 

Given the geometric relationships in Fig. |5} the minimum energy func- 
tion can be written as 

" 2 [(2 + 4 sm\e/2)) /3 + 0.3] ~ ^ ~ 2 sin(0/2) ^^^^ 
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Resting 
Configurations 



Mixed 
Configurations 



Minimum Energy Conjigurations 



Lagrange Resting 

ceo 

Euler Resting 
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Aligned Mixed 
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Transverse Mixed 



Orbiting 
Configurations 



C 
CO 

Lagrange Orbiting 



COO 

Euler Orbiting 



V Resting 



Figure 4: Diagram of all relative configurations for the Full 3-body problem. 
Those boxed can also be stable, and each takes a turn as the global minimum 
energy configuration for a range of angular momentum. 




Figure 5: Generic description of the planar contact geometry between three 
equal sized particles (left) and the mixed configuration geometry (right). 

where the S subscript stands for "Static." The first variation is then 



5£q = sin 9 



3H' 



+ 



(2.9 + 4sin2(0/2))' 8sin^(^/2) 



59 (78) 
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The second variations will be considered on a case- by-case basis. Since 
this system has a constraint on the angle 9^ both the free variations of 
9 and the constrained variations when at the limit must be considered. 



Euler Rest Configuration If ^ = the minimum energy function 
will be stationary. Define this as the Euler Rest Configuration, which 
consists of all three particles lying in a single hne. The stabihty of this 
configuration is evaluated by taking the second variation of the energy 
function and evaluating it at 6' = 180°, yielding 

Recall that the stability condition is that the second variation of the 
energy be positive definite, yielding an explicit condition for stability 
as H'^ > 1.98375, with lower values of H'^ being definitely unstable. 
The Euler Rest Configuration energy can be specified as a function of 
angular momentum: 

H'' _ 5 
2 (2.3) ~ 2 



SE 



(80) 



Lagrange Rest Configuration Now consider the constrained sta- 
tionary point with 9 = 60° (300°). Define this as the Lagrange Rest 
Configuration. Here it suffices to evaluate the first variation at the 
boundary condition, yielding 



SS. 



5 1 60° 



^3 
2 



3H^ 



+ 1 



69 



(81) 



(3.9)2 

At the 60° constraint 69 >0 and the Lagrange Rest Configuration will 
exist and be stable for < 5.07, but beyond this limit an increase in 
9 will lead to a decrease in energy and the relative equilibrium will no 
longer exist. Note that if the 9 = 300° limit is taken, the sign of the 
first variation switches but the constraint surface is now 69 < 0, and 
the same results hold. 

The Lagrange Rest Configuration energy can be specified as a function 
of angular momentum: 



-SL — 



2(1.3) 



(82) 
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Comparing the energy of these two rest configurations shows that the 
Lagrange configuration has lower energy for < 2.99 while the Euler 
configuration has a lower energy above this level of angular momentum. 

2. Variable Contact Configurations 

In addition to the static resting configurations it is also possible to have 
full contact configurations which change as the angular momentum 
varies. These are not fully static, as they depend on having a specific 
level of angular momentum, generating centripetal accelerations that 
balance the gravitational and contact forces. For these configurations, 
as the level of angular momentum varies the configuration itself shifts, 
adjusting to the new environment. For the = 3 case there is only 
one such "variable contact" configuration when restricted to the plane. 
This particular configuration is always unstable, yet plays an important 
role in mediating the stability of the other configurations. 



"V" Rest Configuration The variable contact configuration that 
can exist for this system, called the "V" Configuration for obvious 
reasons, yields the final way for a stationary value of the minimum 
energy function to exist, with the terms within the parenthesis of Eqn. 



78] equaling zero. Instead of solving the resulting quartic equation in 
sin(^/2) it is simpler to evaluate the angular momentum as a function 
of the system configuration to find 

_ (2.9 + 4sin^(g/2))^ 
^ - 24sin=^(^^/2) ^^^^ 

The range of angular momenta that correspond to this configuration 
can be traced out by following the degree of freedom 9 over its range 
of definition. Thus the V Rest Configuration will exist for angular 
momentum values ranging from = 1.98375 at = 180° to = 5.07 
at ^ = 60°. Note that the angles progress from 6 = 180° — )■ 60° as 
the angular momentum increases, and that the limiting values occur 
when the Euler Rest Configuration stabilizes and the Lagrange Rest 
Configuration destabilizes. Note that a symmetric family moves from 
6 = 180° to ^ = 300° at the same levels of angular momentum. 

Taking the second variation and evaluating the sign of along the 
V configuration for arbitrary variations shows that it is always nega- 
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tive definite over the allowable values of 6 and thus that the V Rest 
Configuration is always unstable. 



3. Mixed Configurations 

Now consider mixed configurations where both resting and orbital states 
can co-exist. For = 3 there is only one fundamental topology of this 
class allowed, two particles rest on each other and the third orbits. 
Further, from simple symmetry arguments two candidate states for rel- 
ative equilibrium can be identified, a Transverse Configuration where 
the line joining the two resting particles is orthogonal to the third par- 
ticle {6 = ±90°), and an Aligned Configuration where a single line joins 
all of the mass centers {6 = 0, 180°). To enable a stability analysis a full 
configuration description of these systems is introduced which requires 
two coordinates: the distance from the center of the resting pair to the 
center of the third particle to be R, and the angle between the line R 
and the line joining the resting pair as 9 (see Fig. [s]). 

The distances between the different components can be worked out as 

ri2 = 1 (84) 
= -i? cos + 0.25 (85) 

rgi = + i? cos + 0.25 (86) 

Thus the minimum energy function takes on the form 
" 4(1.2 + i?2) -1 

(87) 



VR^ - RcosO + 0.25 VR^ + RcosO + 0.25 

where the M stands for "Mixed." Taking the variation with respect to 
9 yields 

^ ^ i^sin^ . . 

SeSM = X (88) 
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(i?2 - cos + 0.25)^/^ (i?2 + Rcose + 0.25)^/^ 



As expected, the variation is stationary for the Aligned Configuration, 
6 = 0, 180°, and for the Transverse Configuration, 6 = ±90°. The 
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variation in the distance yields 



M 



3H^R 2R - cos 9 

+ 



2(1.2 + i?2)2 2 (R2 -Rcose + 0.25 f 

2^ + ^°^^ 5R (89) 



2(i?2 + i? cos ^ + 0.25)^/^ 
and is discussed in the following. 



Transverse Configurations First consider the Transverse Config- 
urations with 6 = ±90°. Evaluating the variation of Sm with respect 
to R, setting this to zero, and substituting 9 = 90° allows us to solve 
for the angular momentum explicitly as a function of the separation 
distance 



HI, . ^iH±i?!)L (90) 

3 (i?2 + 0.25)'/' 



This function has a minimum value of angular momentum of Hl^j. ~ 
4.002 . . . which occurs at i? = a/2.6. This is an allowable value of sepa- 
ration and thus this bifurcation will indeed occur. For higher values of 
angular momentum there are two relative equilibria, one with separa- 
tion less than \/2.6 and the other with separation larger than this. The 
inner solution touches the other two particles, forming a Lagrange-like 
configuration, when R = a/3/2. Substituting this into the above equa- 
tion for HIjj, shows that this occurs at a value of 5.07, which is precisely 
the value at which the Lagrange Rest Configuration becomes unstable. 
Recall that this was also the value of angular momentum at which point 
the V Rest Configuration terminated by reaching 60°. Thus at this 
value, which is also equal to the Lagrange Orbit Configuration angular 
momentum at this distance, the inner Transverse Configuration family 
of solutions terminates. Conversely, the outer Transverse Configuration 
persists for all angular momentum values above the bifurcation level. 

Now consider the energetic stability of this class of relative equilibria. 
First note that the cross partials, Sqj^Sm are identically equal to zero 
for the Transverse Configuration. This can be easily seen by taking 



partials of Eqn. 89 with respect to 9 and inserting the nominal value 
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6 = ±90°. Next, taking the second partial of Eqn. 89 with respect to 
6 and evaluating it at the nominal configuration yields 



See£MT = .^o i^df (91) 

and SqqSmt < 0. It is not necessary to check further as this tells 
us that none of the Transverse Configurations are energetically stable. 
The explicit energy of the Transverse Configurations is 



Aligned Configurations Now consider the Aligned Configurations 
with = 0, 180°. Again solve for the angular momentum as a function 
of separation 



ff2 



2(1.2 + 
3^ 



2\2 



1 



(93) 



Finding the minimum point of this equation as a function of R yields 
a cubic equation in R^ without a simple factorization. Root finding 
shows that it bifurcates into existence at a distance of i? = 2.33696 . . . 
with a value of = 5.32417. . .. Again, there is an inner and an 

outer solution. The inner solution continues down to a distance of 
R = 3/2, where the two groups touch and form an Euler configuration. 
The value of the angular momentum at this point equals 6.6125 and 
equals the value at which the Euler Rest Configuration terminates and 
the Euler Orbit Configuration is born. The outer solution continues its 
growth with increasing angular momentum. 

Now consider the energetic stability of these solutions. Similar to the 
Transverse Configurations, the mixed partials of the minimum energy 
function are identically zero at these relative equilibria. The second 
partials of Eqn. |89] with respect to 9 yields 



Sis 



AM 



R 

2" 



(i?- 0.5)3 (i? + o.5)= 



(94) 
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which is always positive. The second partial of Eqn. |89]with respect to 
R is 



SrrSaM 



^ - 0.4) 



4(1.2 + i?2)3 

1 



(i?- 0.5)3 (i? + o.5)= 



{SRf (95) 



The resulting polynomial is of high order and is not analyzed. Alter- 
nately, inspecting the graph of this function shows that it crosses from 
negative to positive at the bifurcation point, as expected. Thus, ini- 
tially the outer Aligned Configurations are energetically stable while 
the inner Aligned Configurations are unstable, and remain so until 
they terminate at the Euler configuration. To make a final check, eval- 
uate the asymptotic sign of the second energy variation. For i? S> 1, 
^MA ~ 4:/ SR. Substituting this into the above and allowing R ^ 1 
again yields dj^^SAM ~ l/R^SR^, and thus the outer relative equilibria 
remain stable from their bifurcation on. 

The explicit energy of the Aligned Configurations are 

c _ 1 1 

" 4(1.2 + i?2) (i?-0.5) " (i? + 0.5) ^^^^ 

A direct comparison between Sma and Smt at the same levels of angular 
momentum shows that the Aligned Configurations always have a lower 
energy than the Transverse Configurations. This is wholly consistent 
with the energetic stability results found throughout. 

4. Purely Orbital Configurations 

Finally consider the purely orbital configurations for this case. As this 
is the sphere restricted problem, the orbital relative equilibria will be 
the same as exist for the point mass problem. 



Euler Solution For the Euler solution take the configuration where 
ri2 = = R and r^i = 2R, R> 1, reducing the configuration to one 
degree of freedom. The minimum energy function then simplifies to 

_ 3i/2 5 

'''' ~ 2(6i?2 + o.9) "2i? ^^^^ 
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Taking the variation of the minimum energy function with respect to 
this configuration then yields 



SrSoe = 9 + w7^ (98) 

Set this equal to zero and solving for the corresponding angular mo- 
mentum to find 

_ 5 i6R^ + 0.9)^ 

It can be shown that there are two orbital Euler configurations for 
H^^ > 8\/5/3 and none for lower values. The non-existence of solutions 
at a given total angular momentum occurs due to the coupling of the 
rotational angular momentum of the different bodies. In this case, 
however, the lower sohitions all exist at i? < 1 and thus are not real for 
this system. In fact, given the constraint R > 1 there will be a single 
family of orbital Euler solutions at Hq^ > 6.6125 with corresponding 
radii ranging from = 1 — )■ oo as Hq^ = 6.6125 — > oo. The correspond 
energy of these Euler solutions as a function of R is 

SoE = -^(6^' -0-9) (100) 

Our simple derivation of the orbital Euler solutions only considers one- 
dimensional variations in the distance. However for a complete stability 
analysis it would be necessary to consider variations of each component 
in turn, as the instability in the Euler problem occurs transverse to the 
line of syzygies. However, the instability of the Euler solutions are well 
known and do not need to be derived here. 



Lagrange Solution To find the conditions for the Lagrange solution 
take the configuration to be ri2 = = r^i = R > 1, again reducing 
the minimum energy function to a single degree of freedom. 

^ 2 i3R^+ 0.9) " I 
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The variation now yields the condition 



3i? 



1 



3H^ 
(3i?2 + o.9)^ 



(102) 



which can be solved for the angular momentum of the orbital Lagrange 
solutions as a function of orbit size 



OL 



3^3 



(103) 



Again, two solutions exist for Hq^ > 16/-\/T0, however the inner solu- 
tion has radius R < 1 and is not allowed by this model. Thus, again for 
the constraint R> 1 there is a single family of Lagrange solution orbits 
that range from i? = 1 — )■ oo as Hqj^ = 5.07 — )■ oo. The corresponding 
energy of these Lagrange solutions as a function of R is 



OL 



1 



{3R^ - 0.9) 



(104) 



Note again that the Lagrange solutions are always energetically unsta- 
ble, and in particular that the equal mass Lagrange solutions are also 
spectrally unstable. As a final point, note that the energy of the Eu- 
ler solutions is actually less than the energy of the Lagrange solutions 
when R^ < 63/60, i.e., when R is near unity. For larger values of R the 
Lagrange solution is always lower energy. 

□ 

Corollary 3. The complete bifurcation chart of relative equilibria, mini- 
mum energy states, and global minimum energy states of the sphere restricted 
N = 3 full body problem as a function of angular momentum follows. It is 
graphically illustrated in Figure\^ 

< < 1.98375 The Lagrange Rest Configuration is the only stable relative 
equilibria and is thus the minimum energy configuration. The Euler 
Rest Configuration does exist, but is unstable. 

= 1.98375 The V Rest Configuration bifurcates into existence from the 
Euler Rest Configuration, and remains unstable throughout its life. The 
Euler Rest Configuration becomes marginally stable. 
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1.98375 < H"^ < 2.99 The Euler Rest Configuration becomes energetically sta- 
ble, but the Lagrange Rest Configuration remains the global minimum 
energy configuration. 

H"^ — 2.99 The Euler and Lagrange Rest Configurations have equal energy 
and are both the global minimum energy configurations. 

2.99 < H'^ < 5.07 The Euler Rest Configuration becomes the global minimum 
energy configuration, while the Lagrange Rest Configuration retains its 
energetic stability. 

H"^ — 4.002 . . . The Transverse Mixed Configurations bifurcate into existence 



at a distance of R — v2.6. 



4.002 . . . < H"^ < 5.07 The Transverse Mixed Configurations continue into 

an inner and outer family at higher values of angular momentum. Both 
of these families are energetically unstable throughout their life. The 
outer family exists for all higher values of angular momentum. 

— 5.07 The V Rest and the Inner Transverse Mixed Configurations ter- 
minate at the Lagrange Rest Configuration. The Lagrange Rest Config- 
uration ceases to exist and transitions into the Lagrange Orbital Con- 
figuration. 

5.07 < < 5.32417. . . The Euler Rest Configuration is the only stable rel- 
ative equilibrium and is thus the global minimum energy configuration. 
The Lagrange Orbital Configuration exists for all higher angular mo- 
mentum and is always unstable. 

= 5.32417. . . The Aligned Mixed Configurations bifurcate into existence 
at a distance of R — 2.33696 .... 

5.32417. . . < H"^ < 5.65907. . . The Inner and Outer Aligned Mixed Config- 
urations exist. The inner configuration is always energetically unstable 
throughout its life. The outer configuration is energetically stable and 
exists for all higher values of angular momentum. The Euler Rest Con- 
figuration remains the global minimum energy configuration. 

H'^ — 5.65907 . . . The Euler Rest and Outer Aligned Mixed Configurations 
have equal energy and are both the global minimum energy configura- 




tions. 
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5.65907. . . < < 6.6125 The Outer Aligned Mixed Configuration becomes 
the global minimum energy configuration, while the Euler Rest Config- 
uration retains its energetic stability. 

= 6.6125 The Inner Aligned Mixed Configurations terminates at the Eu- 
ler Rest Configuration. The Euler Rest Configuration ceases to exist 
and transitions into the Euler Orbital Configuration. 

6.6125 < The Outer Aligned Mixed Configuration is the only stable rela- 
tive equilibrium and is thus the global minimum energy configuration, 
and remains so for all higher angular momentum. The Euler Orbital 
Configuration exists for all higher angular momentum and is always 
unstable. 



1.98375 



4.002.. 



5.07 



o 



I 



o 



• 

cog 



5.07 



5.32417 



6.6125 



Figure 6: Graphic depiction of all relative equilibria and their interactions. 
Green indicates stability while red indicates instability. The global minimum 
is not denoted in this diagram. 



Proof. The transitions can all be evaluated by comparison between the ex- 
istence and stability transitions outlined in the previous proof and by direct 
comparison of energy levels. □ 
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Corollary 4. The N = 3 global and local minimum energy transition dia- 
gram as a function of angular momentum is discontinuous and exhibits hys- 
teresis. Figure \2\ tracks the energetically stable configurations as a function 
of angular momentum, shows where explicit transitions occur, and indicates 
the excess energy at each transition. To track the evolution of a system un- 
der changing angular momentum will require dynamical analysis at specific 
transition points. 




Figure 7: Bifurcation diagram showing the energy-angular momentum curves 
of all stable relative equilibria and their transition paths for increasing and 
decreasing angular momentum. 



Proof. The plot is created by comparing the energies of the different energet- 
ically stable configurations and computing energy differences at transitions. 
Since excess energies exist at different transition points, it will be necessary 
for the dynamical evolution of the system to be followed to determine actual 
outcomes. □ 



4.4 iV > 1 

The paper ends with a few observations for when is a large number. Sim- 
ulations of large numbers of grains that gravitationally attract each other 
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and can rest on each other have been explored over the last decade using a 
variety of simulation methods. The original work in this area was performed 
by Richardson and collaborators (c.f. and used hard-sphere models for 
surface interactions. More recent models using soft-sphere models for grain 
interactions have been explored by Sanchez and Scheeres [T3], and in some 
cases have advantages in terms of tracking the system-wide energy, angular 
momentum, and internal forces of these systems when there are significant 
components that rest on each other. Investigations of relative equilibria of 
the minimum energy function for such large particle simulations could be 
investigated, but have not been as of yet. 

There are several direct observations that can be made for large sys- 
tems. First, for non-rotating systems it is expected that the global minimum 
configuration will tend towards a spherical collection of particles resting on 
each other. Depending on the presence and strength of friction between 
the grains, it is also possible that a variety of non-spherical shapes may 
be stable relative equilibria, thus generalizing the ideal results from contin- 
uum mechanics and theory. Likewise, as angular momentum is increased the 
classical Maclaurin and Jacobi ellipsoids should always be relative equilib- 
ria, however one must again note the probable existence of multiple other 
non-ellipsoidal configurations that are also relative equilibria. This is to be 
expected, given the presence of such multiple configurations at a given level 
of angular momentum for particle systems with N > 2. Of specific inter- 
est are the evolutionary pathways and transitions for these systems as the 
total angular momentum is increased (or decreased) mimicking the natural 
evolution of rubble pile asteroids under the YORP effect. Analytical studies 
of this have been made by Holsapple |H] and recreated using simulations by 
Sanchez [TB] . Of particular interest for these studies are the pathways to fis- 
sion of large collections of particles. A systematic study of these transitions 
will be of real interest in understanding the overall evolution and creation of 
binary bodies among small asteroids, at the least [9]. 

To close this section, a final hypothesis is given. 

Hypothesis 1. For all N > 2, energetically stable relative equilibrium con- 
figurations of Full Body problems will result in one of two general states: 

1. All the particles rest on each other in one collection and spin at a 
uniform rate; 

2. All the particles separate into two collections in a mutually circular 
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orbit about each other with doubly-synchronous rotation. 

Further, that all configurations where the particles condensate into 3 or more 
collections will be energetically unstable and will either eject the excess collec- 
tions or, under energy dissipation, will condense into one or two collections. 

One possible avenue to the proof of this hypothesis is through the proof for 
the point mass N body problem that all orbital relative equilibria with N >3 
are energetically unstable. As for most collections with > 3 in relative 
equilibria, the leading order of their gravitational attractions will always 
represent the point mass body problem, with higher-order gravitational 
perturbations arising from the non-spherical configuration of the condensed 
grains in each collection. 

5 Conclusions 

This paper considers the basic and fundamental question: What are the min- 
imum energy and stable configurations of the A^-body problem in celestial 
mechanics for a given value of angular momentum. It is shown that this 
problem is ill-posed for the classical point mass model of celestial mechanics, 
but is well posed if all density distributions for these systems are made finite. 
This stipulation requires that rotational energy and angular momentum also 
be considered in searching for these stable configurations, and indeed funda- 
mentally changes the problem. Specifically, under this reasonable stipulation 
the minimum energy solutions of the finite density 2-body problem change 
fundamentally as compared to the point mass 2-body problem. This ques- 
tion is also rigorously explored for A^ = 3 along with a number of hypotheses 
for the A^ ^ 1 problem. In the course of these discussions a number of 
fundamental questions are identified on the evolution of celestial mechanics 
systems under an increase of angular momentum, as has been documented 
to occur in nature. 
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